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Figure 1: A sketch of the QCD Phase Diagram, according to the collaboration CBM at FAIR [[l|], it is 
usually assumed that the critical point for deconfinement coincides with the critical point for chiral symmetry 
restoration. 



1. Introduction 

Our main motivation is to contribute to understand the QCD phase diagram for finite T and 
/I. The QCD phase diagram is scheduled to be studied at LHC, RHIC and FAIR, and is sketched 
in Fig |l|. Notice that, after enormous theoretical efforts, the analytic crossover nature of the finite- 
temperature QCD transition was finally determined by Y. Aoki et al. ||^, utilizing Lattice QCD 
and physical quark, and reaching the continuum extrapolation with a finite volume analysis. 

Here we utilize the Coulomb gauge hamiltonian formalism of QCD, presently the only con- 
tinuum model of QCD able to microscopically include both a quark-antiquark confining potential 
and a vacuum condensate of quark-antiquark pairs. This model is able to address excited hadrons 
as in Fig. ^ and chiral symmetry at the same token, and we recently suggested that the infrared 
enhancement of the quark mass can be observed in the excited baryon spectrum at CBELSA and 
at JLAB |j3|, 1^. Thus the present work, not only addresses the QCD phase diagram, but it also 
constitutes the first step to allow us in the future to extend the computation of nay hadron spectrum, 
say the Fig. ^computed in reference |^, to finite T . 

Both for the study of the hadron spectra, and for the study the QCD phase diagram, a finite 
quark mass is relevant. In the phase diagram, a finite current quark mass mo affects the position 
of the critical point between the crossover at low chemical potential pL and the phase transition at 
higher /i. Moreover the current quark mass affects the QCD vacuum energy density relevant 
for the dark energy of cosmology. This all occurs in the dynamical generation of the quark mass 
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Figure 2: First calculation of excited baryons with a chiral invariant quark model | 



m{p). While the quark condensate (yy) is a frequently used order parameter for chiral symmetry 
breaking, the mass gap, /. e. the quark mass at vanishing momentum m(0) is another possible 
order parameter for chiral symmetry breaking. 

Here we address the finite temperature string tension, the quark mass gap for a finite current 
quark mass and temperature, and the deconfinement and chiral restoration crossovers. We conclude 
on the separation of the critical point for chiral symmetry restoration from the critical point for 
deconfinement. 



2. Fits for the finite T string tension from the Lattice QCD energy Fi 

At vanishing temperature T = 0, the confinement, due to the formation of chromo-electric 
and chromo-magnetic flux tubes as shown in Figs. || and ^ can be modelled by a string tension, 
dominant at moderate distances, 

V{r)^^ + Vo + ar. (2.1) 

At short distances we have the Luscher or Nambu-Gotto Coulomb potential due to the string vi- 
bration plus the One Gluon Exchange Coulomb potential, however the Coulomb potential is not 
important for chiral symmetry breaking. At finite temperature the string tension g{T) should also 
dominate chiral symmetry breaking, and thus one of our crucial steps here is the fit of the string ten- 
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Figure 3: The fields generated by static hybrid sources, notice how the flux organizes into fundamental 
colour triplet flux tubes. 

sion <j{T) obtained from the Lattice QCD data for the quark-antiquark free energy of the Bielefeld 
Lattice QCD group, [|,| ^ I g]. 

At finite temperature, the quark, or the quark-antiquark free energies, can be computed utiliz- 
ing Polyakov loops, as in Fig. |[ The Polyakov loop is a gluonic path, closed in the imaginary time 
t4 (proportional to the inverse temperature T^^) direction in a periodic boundary Euclidian Lattice 
discretization of QCD. They measures the free energy F of one or more static quarks. 



P{0) =Ne 



P"{0)P"{r) =Ne 



(2.2) 



If we consider a single solitary quark in the universe, in the confining phase, his string will travel as 
far as needed to connect the quark to an antiquark, resulting in an infinite energy F. Thus the 1 quark 
Polyakov loop f is a frequently used order parameter for deconfinement. With the string tension 
a(r) extracted from the qq pair of Polyakov loops we can also estimate the 1 quark Polyakov loop 
P(0). 

At finite T , we use as thermodynamic potentials the free energy F\ , computed in Lattice QCD 
with the Polyakov loops [|5|, ^ 0, ||, and illustrated in Fig. ^. It is related to the static potential 
V{r) = —fdr with F\{r) = —fdr — SdT adequate for isothermic transformations. In Fig. ^we 
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Figure 4: The perpendicular profile of the longitudinal electric field, in units of the lattice spacing a — 
0.07261 (85)/m. The flux tube is so thin, that in quark models it is usually modelled by a single parameter, 
the string tension a. 



extract the string tensions g{T) from the free energy Fi (T) computed by the Bielefeld group, and 
we also include string tensions previously computed by the Bielefeld group [[Io|]. 

We also find an ansatz for the string tension curve, among the order parameter curves of other 
physical systems related to confinement, i. e. in ferromagnetic materials, in the Ising model, in 
superconductors either in the BCS model or in the Ginzburg-Landau model, or in string models, 
to suggest ansatze for the string tension curve. We find that the order parameter curve that best 



fits our string tension curve is the spontaneous magnetization of a ferromagnet [11], solution of the 
algebraic equation, 



M 



tanh 



Tc M 

TmZ, 



(2.3) 



In Fig. ^ we show the solution of Eq. g3| obtained with the fixed point expansion, and compare it 
with the string tensions computed from lattice QCD data. 



3. The mass gap equation with finite T and finite current quark mass mo. 



Now, the critical point occurs when the phase transition changes to a crossover, and the 
crossover in QCD is produced by the finite current quark mass mO, since it affects the order pa- 
rameters P or a, and the mass gap m(0) or the quark condensate {qq). Moreover utilizing as order 
parameter the mass gap, i. e. the quark mass at vanishing moment, a finite quark mass transforms 
the chiral symmetry breaking from a phase transition into a crossover. For the study of the QCD 
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Figure 5: The renormalized Polyakov loop of a single quark, computed by our Lattice QCD group. 

phase diagram it thus is relevant to determine how the current quark mass affects chiral symme- 
try breaking, in particular we study in detail the effect of the finite current quark mass on chiral 
symmetry breaking, in the framework of truncated Coulomb gauge QCD with a linear confining 
quark-antiquark potential. 

The most fundamental information for the quark-antiquark interaction in QCD comes from 
the Wilson loop in Lattice QCD, providing the confining quark-antiquark potential for a static 
quark-antiquark pair. This potential is consistent with the funnel potential, also utilized in the 
quark model to describe the quark-antiquark sector of meson spectrum, in particular to describe 
the linear behaviour of mesonic Regge trajectories. Notice that the short range Coulomb potential 
could also be included in the interaction, but here we ignore it since it only affects the quark mass 
through ultraviolet renormalization [|2|], which is assumed to be already included in the current 
quark mass. Here we specialize in computing different aspects of chiral symmetry breaking with 
linear confinement V = or. Since we are interested in working at finite temperature T we utilize a 
recent fit of lattice QCQ data with a temperature dependent string tension o{T). 

To address the light quark sector it is not sufficient to know the static quark-antiquark potential, 
we also need to know what Dirac vertex to use in the quark-antiquark-interaction. This vertex is 
necessary to study not only the meson spectrum but also the dynamical spontaneous breaking of 
chiral symmetry. To determine what vertex to use, we review how the quark-antiquark potential 



can be approximately derived from QCD, in two different gauges. In Coulomb gauge [|130, 

V-A(x,0 =0 



(3.1) 



the interaction potential, as derived by Szczepaniak and Swanson [14, [T^], is a density-density 
interaction, with Dirac structure (8 y". Another approximate path from QCD considers the mod- 
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Figure 6: We show examples finite temperature static quark-antiquark potentials, in particular the T < Tc 
and T > Tc Lattice QCD data for the free energy Fi, thanks to ^ ^ |[ ^ Olaf Kaczmarek et al. The solid 
line represents the T = static quark-antiquark potential. In this paper we discuss the use of the free energy 
as a finite temperature quark-antiquark potential. 
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Figure 7: Detail of the string tension fit in the case of J = 0.947^. We cut the low distance part in such a 
way that a linear fit is stable for cuttof changes. 



ified coordinate gauge of Balitsky Q16| ] and in the interaction potential for the quark sector, retains 
the first cumulant order, of two gluons 



17|, |18), |19||. This again results in a simple density-density 
effective 7*^ (g) /* confining interaction. As in QCD, this only has one scale, say Kq, in the inter- 
action, since both the quark condensate and the hadron spectrum turn out to be insensitive to any 
constant U in the potential. Thus our framework is similar to an expansion of the QCD interaction, 
truncated to the leading density-density term, where the confining quark-antiquark potential is a 
linear potential. 

While this is not exactly equivalent to QCD, our framework maintains three interesting aspects 
of non-pertubative QCD, a chiral invariant quark-antiquark interaction, the cancellation of infrared 
divergences [20,2l|, 23, 24, 25], and a quark-antiquark linear potential [26, 27, | 



14, 2 



Importantly, since our model is well defined and solvable, it can be used as a simpler model than 
QCD, and yet qualitatively correct, to address different aspects of hadronic physics. In particular 
here we study how chiral symmetry breaking occurs at finite temperature T and chemical potential 
jj., in the realistic case of small but finite current quark masses. Thus we apply our framework to 
the phase diagram of QCD. 

Our interaction potential for the quark sector is. 



Vi 



{mop 
VAnx)A>(x)^V( 



- ia ■ V) (//-(x) + ^ y d'^y 
|x-y|)i/At(y)AV(y) 



(3.2) 



where the density-density interaction includes just the linear confining potential together with an 
infrared constant, which may be possibly divergent. 

The mass gap equation and the energy of a quark are determined from the Schwinger-Dyson 
equation at one loop order using the Hamiltonian of Eq. (3.2). for a recent derivation with all 
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Figure 8: (top) The critical curve for -^-^ as a function of |r, for T ~ Tc it behaves like a square root, 
(bottom) Comparing the magnetization critical curve with the string tension (j/(Jq, fitted from the long 
distance part of Fi , they are quite close. 



details see [plQ. The interaction in the four momentum of the potential and quark propagator term 
includes an integral in the energy 



1 



2n p°-E{p) + ie 2 ' 



(3.3) 
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Figure 9: The linear potential must be regularized for the Fourier Transform, and we show our two different 
regularizations both leading to V{r) — > when r — > 0. Two successions of curves are plotted, leading exactly 
to linear potentials in the limit of a vanishing infrared regulator ji The regularization of the negative 
curves maintains the potential monotonous, but adds an infrared negative constant to the potential. The 
regularization of the positive curves maintains V (0) =0 but the potential decreases for large r. 



which factorizes trivially from the vector p momentum integral. Using spherical coordinates, the 
angular integrals can be performed analytically and finally only an integral in the modulus of the 
momentum remains to be computed numerically. We arrive at the mass gap equation in two equiv- 
alent forms, of a non-linear integral functional equation, 

(T f°° flk 

= pS{p)-m,C{p)-^ / ^[ (3.4) 

Jo 2n '- 

lA{k,p,^)S{k)C{p)-lB{k,p,^)S{p)C{k)] , 
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and of a minimum equation of the energy density (o, 

2p^C{p) + 2p^moS{p) + a X 
Jo ^^A{k,p,^)S{k)S{p) +lB{k,p,^)C{p)C{k) 



In Jo 2k 
dk 



(3.5) 



where the functions /g and Ia are angular integrals of the Fourier transform of the potential. In what 
concerns the one quark energy we get, 



E{p) = pC{p)+moS{p) + - 



dk 

p^ Jo 2n 



lA{k,p,}i) X 



S{k)S{p)+lB{k,p,pL)C{p)C{k) 



(3.6) 



In the chiral limit of massless current quarks, the breaking of chiral symmetry is spontaneous. 
But for a finite current quark mass, some dynamical symmetry breaking continues to add to the 
explicit breaking caused by the quark mass. The mass gap equation at the ladder/rainbow truncation 
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Figure 11: We plot the mass gap m(0) as a function of T for the hght u (hghtest) and d (shghtly heavier) 
quarks. We have a crossover, close to a phase transition, since the current masses of the light quarks are 
much smaller than the dynamically generated constituent quark mass. 



of Coulomb Gauge QCD in equal time reads, 



m{p) = mo + 
h{p,k,ix) = 
lB{p,k,n) = 



dk lA{p,k,lx)m{k)p — lB{p,k,lJ.)m{p)k 



p^ Jo 271 
pk 



y/k^ + m{kY 
pk 



{p-kf + jx^ [p + kf + jx'^ 

pk pk 1 {p-kf + jx^ 

{p - kY +ix^^ {p+ky+^^^2 °^ (p + k)^ + 



(3.7) 



The mass gap equation ( 3.7 ) for the running mass m{p) is a non-linear integral equation with a nasty 



cancellation of Infrared divergences [pA 33, 34 1. We devise a new method with a rational ansatz, 



and with relaxation [31 1, to get a maximum precision in the IR where the equation extremely large 
cancellations occur. Since the current quark masses of the six standard flavours u, d, s, c, b, t span 
over five orders of magnitude from 1.5 MeV to 17 1 GeV, we develop an accurate numerical method 
to study the running quark mass gap and the quark vacuum energy density from very small to very 
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Figure 12: Log plot of the mass gap mc/{0) for the six different flavours of quarks as a function of T. From 
bottom to top we show the quarks u, d, s, c, b and f . The heavier the quark, the weaker the crossover gets. 



large current quark masses. The solution m{p) is shown in Fig. [T^ for a vanishing momentum 
p = 0. 

At finite T , one only has to change the string tension to the finite T string tension a(r) of Eq. 



{13), for different quark masses [35|, and also to replace an integral in by a discrete sum in 
Matsubara Frequencies. Both are equivalent to a reduction in the string tension, a — a* and thus 
all we have to do is to solve the mass gap equation in units of o* . The results are depicted in Fig. 
10 . Thus at vanishing mo we have a chiral symmetry phase transition, and at finite wiq we have a 
crossover, that gets weaker and weaker when mq increases. This is also sketched in Fig. 



4. Infrared regularization of the linear confining potential and the Matsubara sum 

Notice that in the case of a linear potential, divergent in the infrared, the Fourier transform 
needs an infrared regulator jx eventually vanishing. We illustrate two possible infrared regulariza- 
tions of the linear potential if Fig. ^ 

A possible regularization of the linear potential is. 



V{r) 



At 



~ 1- or , 

At 



(4.1) 



corresponding to a model of confinement where the quark-anti quark system has an infinite binding 
energy ^at the origin r = 0, is monotonous and only vanishes at an arbitrarily large distance. This 



13 



quark mass and meson spectrum at deconfinement 



Pedro Bicudo 



FrCr) 
2m( 



nip » c 



2 m, 



iHq finite 



saturation 
at T<Tc 



2mo 



mo=0 



Figure 13: Sketch of the saturation of confinement at the energy of two heavy-light groundstate mesons. 

potential has a simple three-dimensional Fourier transform, 

47rrsin(^r) 



Jo k 



(4^2, 



and this is the most common form of the linear potential in momentum space utilized in the lit- 
erature. Notice that this is infrared divergent due to the — infinite binding energy in the limit 
where the regulator jx ^0. 

If we want to avoid the infinite binding energy we should use a different regularization of the 
linear potential, also vanishing when r ^ oo but not monotonous since it grows linearly at the origin 
starting with V(0) =0, 

V{r) = are-^'' (4.3) 

where the Fourier transform, 

is such that the integrals in k no longer diverge. For instance, 

/ + 00 
k^dkV{k)=0 (4.5) 
-oo 
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Figure 14: Sketch the order parameter P polyakov loop for the quark antiquark-system, as impUed by the 
saturation of the energy of the quark-antiquark system. 



since this is proportional to V{0) = 0. The new term in the potential ^li^^i-^^ is equal to {2Ti)^5^{k) / pL 
in the limit jU — 0, and this potential is infrared finite. Both the potentials in Eqs. ( p~I| ) and ( p3| ) 
are illustrated in Fig. ^ 

In the vanishing temperature limit T = the different regularizations lead to the same physical 
results since any constant term in a density-density interaction has no effect in the quark running 
mass m{p) or in the hadron spectrum [^]. The regularizations only contribute to the potential 
and the one-quark energy, but it occurs that these contributions exactly cancel in the chiral order 
parameters, and in the hadron spectrum. 

However at T 7^ the two different regularizations may lead to different physical results. 

In the r = mass gap equation or Schwinger Dyson equation, we have the Minkowski space 
integral in p^ of the quark propagator pole of Eq. (3.3 ) and this is equivalent to an integral in p'^ in 



Euclidian space after a Wick rotation in the Argand space, 

<'+'"idp'^ i 1 



271 ip"^ -E{ip) + ie 2 ■ '■'^■^^ 
real axis the path corresponds to z = p^ and in the imaginary axis the path corresponds to z = ip^- 



Notice that the integral in Eq. is only identical to the one in Eq. ( p.3[ ) if the one quark 
energy E{p) > 0. If E{p) < the integral of Eq. ( pT^ changes sign, and this is consistent with the 
pole moving from the fourth quadrant to the third quadrant of the Argand plane. While the integral 
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T 



Figure 15: We sketch the hght hadron spectrum. Notice that with T ^0, the light masses and wavefunctions 
do essentially scale with <y{T), and thus also essentially vanish at T = . 



in Eq. (3.3 ) is insensitive to this translation of the pole, the Wick rotation leads to a pole correction. 



For simplicity, we choose to work with positive one quark energies only E(p) > 0. 

In finite temperature T and density p, the continuous euclidian space integration of eq. (16) 
is extended to the sum in Matsubara frequencies, 

J^J''^i{2n + \)KkT-[E{p)] ■ ^"^-"^^ 
It is clear that in the vanishing temperature and density kT «E limit one gets back the initial Eu- 



clidean space integral of eq. (16), when the Matsubara sum approaches the continuum integration 



with kT 

Notice that at T 7^ the results depend on the constant in the potential. Since the Matsubara 
sum depends on the one quark energy E{p), and a constant shift of the potential added to the linear 
term, say V = a{T)r—UoUo > 0, affects the energy hy a ^Uq, thus weakening the finite T effects 
scaling like On the other hand if the shift is small enough to produce at some p a negative 
energy E{p) < 0, at vanishing density jJ. =0, then the Matsubara sum could break down. 

Thus we study here only one scenario for the constant potential shift —Uq, present anyway 
in the different possible infrared regularizations of the linear potential. Our scenario consists in a 
maximum energy, with ?7o — > 00 as in the standard regularization of the linear potential of Eq. (p~l|). 
In that case E{p) 00 and the Matsubara sum is simply constant, producing always the T = 



result of "2^. Our results are shown in Figs. 1 1 and 12 



We don't explore this second possible scenario here, leaving it for a future study, consisting in 



having the minimal Uq, closer to the regularization of the linear potential of Eq. ( 4.3 ), just sufficient 



to cancel the one quark energy at vanishing momentum E{0) = 0. In that case the Matsubara 
sum makes a larger difference, interpolating between at vanishing momentum and — 5 at large 
momentum. But we expect that, whatever the contribution of the Matsubara sum is, the temperature 
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Figure 16: We sketch the behaviour of the meson spectra with temperature. The light hadron masses M^^ 
are dominated by the hnear confining potential, and essentially vanish at Tc. On the other hand, the heavy 
hadron masses Mqq are also significantly bound by the Coulomb potential, providing enough binding up to 
a melting temperature above Tc. 



dependence of the string tension (j{T) is the dominant finite T effect, and thus for a first exploration 
our first scenario is sufficient. 



5. Chiral symmetry and confinement crossovers with a finite current quark mass 

We now study whether the two main phase transitions in the QCD phase diagram, confinement 
and chiral symmetry breaking, have two different two critical points or a coincident one. Confine- 
ment drives cfiiral symmetry breaking, and at small density both transitions are a crossover, and 
not a first or second order phase transition due to the finite quark mass. 

However the quark mass affects differently these two phase transitions in the QCD. In what 
concerns confinement, the linear confining quark-antiquark potential saturates when the string 
breaks at the threshold for the creation of a quark-antiquark pair. Thus the free energy F{0) of 
a single static quark is not infinite, it is the energy of the string saturation. The saturation energy is 
of the order of the mass of a meson i. e. of 2mo. For the Polyakov loop we get, 

P(0) -A^e^^"""/^ . (5.1) 

Thus at infinite mo we have a confining phase transition, while at finite mo we have a crossover. 



that gets weaker and weaker when mo decreases. This is sketched in Figs. |13|and |14. 

Since the finite current quark mass affects in opposite ways the crossover for confinement and 
the one for chiral symmetry, we conjecture that at finite T and jJ. there are not only one but two 
critical points (a point where a crossover separates from a phase transition). Since for the light u 
and d quarks the current mass mo is small, we expect the crossover for chiral symmetry restoration 
critical to be closer to the jJ. =0 vertical axis, and the crossover for deconfinement to go deeper 
into the finite n region of the critical curve in the QCD phase diagram depicted in Fig. [T[ 

We now compute the hadron spectrum, in particular the meson spectra. In what concerns the 
light meson masses Mqq are dominated by the linear confining potential. At Tc, the string tension 
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a(r) vanishes, the confining potential disappears, and thus all light hadrons decrease their mass 



until they melt at T = T^.. Notice that chiral symmetry is not entirely restored because mo/y (y{T) 
increases with T. Nevertheless, with T ^ 0, the masses and wavefunctions do essentially scale 



with a(r), and thus also essentially vanish at T = T^. This is sketched in Fig. |T5[ A comparable 
vanishing of light meson masses also occurs in the sigma model [^6|]. 

In what concerns the heavy meson masses Mqq, the Coulomb potential is also relevant, and 
thus the groundstates of the heavy mesons still bind above Tc. Thus while the light mesons melt at 



T = Tc, the heavy meson groundstates melt at a T > [ ]37| , pq], as depicted in Fig. |16 
6. Outlook 

We remark that the pure gauge string tension cj(r) is well fitted by the condensed matter 
physics magnetization curve x, and we utilize it. 

We compute the dynamically generated quark mass m{p) , solving the mass gap equation both 
for finite current quark masses mo and for finite T. The finite current quark masses turn both the 
confinement and the chiral symmetry phase transitions into two different crossovers. 

We qualitatively study the full spectra of light hadrons at finite T, including the excited spectra, 
and conclude that the light hadron masses essentially vanish at T = Tc. The light hadrons all melt 
at r = Tc, since the masses and wavefunctions essentially scale with o{T). 

We soon plan to complete the part of this work which is only sketched here, /. e. to compute the 
excited hadron spectrum at finite T, continuing the work initiated with Tim Van Cauteren, Marco 
Cardoso, Nuno Cardoso and Felipe Llanes-Estrada [^, and to compute and compare the crossover 
curves for the chiral symmetry restoration and for the deconfinement at finite T. 
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